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Abstract 

We analyze the effect of the non-universal Z boson in the rare decays B s —> l + l ~, B s —> HU'y 
and Z —> bs decays. These decays involve the FCNC mediated b —> s transitions, and are found 
to be very small in the standard model. The smallness of these decays in the standard model 
makes them sensitive probe for new physics. We find an enhancement of at least an order in these 
branching ratios because of the non-universal Zbs coupling. 
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The rare decays induced by the flavor changing neutral current (FCNC) transitions are 
very important to probe the flavor sector of the standard model (SM). In the SM they 
arise from one-loop diagrams and hence are generally suppressed in comparison to the tree 
diagrams. This makes them as a sensitive probe for new physics. The rare B decays which 
are mediated by the FCNC transitions are of the kind b —> s or b —> d. Prominent examples 
of rare B decays are B —> K* 7, B —> pq, B —> B S) d —> l + l~. During the last few 

years, considerable theoretical attention has therefore been focused on these decays in view 
of the planned experiments at B-factories, which are likely to measure branching fractions as 


low as 10 8 


AS- 


The results of the branching ratios of FCNC mediated H-decays are very 
sensitive to new physics beyond the SM. Thus, a detailed investigation of the rare decays is 
a promising way to discover or severely constrain the new physics. 

In this paper, we focus on a specific class of rare decay modes B s —> l + l~ and B s —> / + / _ y, 
which are mediated by the Z boson exchange and the rare Z decay mode Z —» bs. We 
consider the effect of the non-universal Z boson which induces FCNC interaction at the tree 
level. It is well known that FCNC coupling of the Z boson can be generated at the tree level 
in various exotic scenarios. Two popular examples discussed in the literature are the models 
with an extra U( 1) symmetry [3] and those with the addition of non-sequential generation 
of quarks 4J]. In the case of extra 17(1) symmetry the FCNC couplings of the Z boson are 
induced by Z — Z' mixing, provided the SM quarks have family non-universal charges under 
the new 17(1) group. In the second case, adding a different number of up- and down-type 
quarks, the pseudo CKM matrix needed to diagonalize the charged currents is no longer 
unitary and this leads to tree level FCNC couplings. It should be noted that, recently, there 
has been renewed interests shown in the literature concering the non-universal Z induced 


new physics |5[. In light of this it necessitates also a detailed investigation of rare B decays, 
which are very promising to discover and/or to constrain new physics. 


Here we will follow the second approach [ 4 ] to analyze some FCNC induced rare decays. 
It is a simple model beyond the standard model with an enlarged matter sector due to an 
additional vector like down quark D±. The presence of an additional down quark implies a 
4x4 matrix V ia (i = u, c, t, 4, a = d, s, 6, b r ), diagonalizing the down quark mass matrix. For 
our purpose the relevant information for the low energy physics is encoded in the extended 
mixing matrix. The charged currents are unchanged except that the Vqkm is now the 3x4 
upper sub-matrix of V. However, the distinctive feature of this model is that the FCNC 
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interaction enters neutral current Lagrangian of the left handed down quarks as 


— 


9 


2 cos 9w 


ULil^u Li - d La U a[ rf ,l d L g - 2 sin 2 6 w J£ n , 


i 


( 1 ) 


with 

Uae = Y. H V ‘D = s «e - v £> v w . ( 2 ) 

i=u,c,t 

where U is the neutral current mixing matrix for the down sector, which is given above. As 
V is not unitary, U ^ 1. In particular the non-diagonal elements do not vanish. 


u a0 = -V 4 * a V 4i3 ^ 0 for a^/3. 


( 3 ) 


Since the various U a p are non vanishing, they would signal new physics and the presence 
of FCNC at the tree level and this can substantially modify the predictions of SM for the 
FCNC processes. 

Now let us consider the FCNC process B s —> l + l~. These decays, in particular the process 
B s —■> /i + p~ has attracted a lot of attention recently since it is very sensitive to the structure 
of SM and potential source of new physics beyond the SM. Furthermore, this process is very 
clean and the only nonperturbative quantity involved is the decay constant of B s meson 
which can be reliably calculated by the well known non-perturbative methods such as QCD 
sum rules, lattice gauge theory etc. Therefore, it provides a good hunting ground to probe 
for new physics. The branching ratio for B s —> l + l~ has been calculated in the SM [a) and 


also in beyond the SM in a number of papers [?]. Let us start by recalling the result for 
B s —> l + l~ in QCD-improved standard model. The effective Hamiltonian describing this 
process is 


Gpoi 

Heff = -4rV tb v u 


V2' 


7T 


c e 9 ff (s 7/t P L b)(l Y l) + C'io (s 7 /, Pl b)(l Y 75 l) 


^l^(sia^P R b)(Jr l) 


( 4 ) 


where Pl,r = | (1 =F 75 ) and q is the momentum transfer. CV^are the Wilson coefficients 
evaluated at the b quark mass scale in NLL order with values 


C e 7 ff = -0.308 , C 9 = 4.154 , C'io = -4.261 . 


( 5 ) 


The coefficient C 9 ^ has a perturbative part and a resonance part which comes from the 
long distance effects due to the conversion of the real cc into the lepton pair l + l~. Hence, 
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Cg^ can be written as 


C e 9 ff = C 9 + Y(s) + C r 9 es , (6) 

where the function Y (s) denotes the perturbative part coming from one loop matrix elements 
of the four quark operators and is given in Ref. j^]. The long distance resonance effect is 
given as [u]] 

m Vi r(Vi -> m~) 


Qyj- 

Cg E " = —— ( 3 C*i + C*2 + 3C3 + C4 + 3C5 + (%) 'y 

cr 


k- 


rnY — s — imv T y. ’ 
Vi 1 1 


(7) 


where the phenomenological parameter k is taken to be 2.3, so as to reproduce the correct 
branching ratio B{B -»• J/ipK* -> /TZ+Z") = B{B -»• J / ipK*)B{J/ip -> Z+Z"). In this 
analysis, we will consider only the contributions arising from two dominant resonances i.e., 
J/ip and ip'. The values of the coefficients Q’s in NLL order are given in 8 ] as Ci = 
-0.151 , C 2 = 1.059 , C 3 = 0.012 , C 4 = -0.034 , C 5 = 0.010 and C 6 = -0.040. 

To evaluate the transition amplitude one can generally adopt the vacuum insertion 
method, where the form factors of the various currents are defined as follows 


(0 | s 7 m 75 b | B° s ) = if Bs PB > 

(0 | s 75 b\B° s ) = if Ba m Bs , 

(0| s P R b\B° s ) = 0 . ( 8 ) 


Since p^ = + P-, the contribution from Cg term in Eq. (JJJ) will vanish upon contraction 

with the lepton bilinear, C 7 will also give zero by (JBJ) and the remaining C 10 term will get a 
factor of 2 m/. 

Thus the transition amplitude for the process is given as 

M(B S - Z+Z") = V tb V* f Ba Cio rm {Hi) , (9) 

v 27T 

and the corresponding branching ratio is given as 


Bin, 


;+ n = Gl r B , 


167T 3 


« 2 f Bs m Ba m\} \V tb V*\- C 


Amf 


10 


m 


( 10 ) 


Ba 


Helicity suppression is reflected by the presence of mf in (Hhd which gives almost vanishingly 
small value for e + e _ and a very small branching ratio of (3.4 ± 0.5) x 10 -9 for /i + U|- 
The published Tevatron/CDF physics results with luminosity 171 pb -1 provides the bound 
on B s —¥ [1 


B{B S fx+fj,-) < 5.8 x 10 ~ 7 (90% C.L.) . 


( 11 ) 
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FIG. 1: Feynman diagram for B s —► l + l in the model with tree level FCNC transitions, where 
the blob represents the tree level flavor changing vertex. 

Recently, this branching ratio has been further constrained by the DO collaboration Q] with 
an upper bound 

B(B S -> < 5.0 x 1(T 7 (95% C.L.) . (12) 

It should be noted that the r channel is free from this hclicity suppression however, its 
experimental detection is quite hard due to the low detection efficiency and that is why we 
do not have any experimental upper limit for this process as yet. 


Now let us analyze the decay modes B s —> l + l in the model with the Z mediated FCNC 


occurring at the tree level 


Q. 


The corresponding diagram is shown in Figure-1, where the 
blob represents the tree level Zbs coupling. The effective Hamiltonian for B s —> l + l~ is 
given as 

Gf tt r-_ u/ ' 1 - „l_ . \ 7 l (13) 


Heff = Usb ~~ 75 ^ - C l jd^)i\ , 

where C l v and C l A are the vector and axial vector Zl + l~ couplings, which are given as 


C l v — — — + 2 sin 2 9 


w 


n l — _ ■ 
u a — ; 


(14) 


Since, the structure of the effective Hamiltonian m in this model is the same form as that 
of the SM, like ~ {V — A)(V — A) form, therefore its effect on the various decay observables 
can be encoded by replacing the SM Wilson coeiffients (CqG^ sai and (Cio) SAI by 


(jeff _ ^rieff^SM _|_ 27T U s bCy ^e// _ _ 2?r UsbC l A 


(15) 


« V tb v t * s 7 iU v a V tb V t * s 

It should be noted that U sb is in general complex and hence it induces the weak phase 
difference (9) between the SM and new physics contributions. Since the value of the Wilson 
coefficients C 9 and C A q are opposite to each other as seen from Eq. ©, and the new 
physics contributions to Cg and C\g are opposite to each other, one will get constructive or 


5 





destructive interference of SM and NP amplitudes for 0 = tt or zero (where 0 denotes the 
relative weak phase between SM and NP contribution in the above equation). 

Thus, one can obtain the branching ratio including the NP contributions by substituting 
from (gSD in (Hill) . Now using the value of \Ui )S \ ~ 10 3 Q], which has been extracted 
from the recent data on B(B —> X s l + l~), sin 2 9 W = 0.23, the particle masses from 


0 . 


a = 1/127, the decay constant f b s = 0.24 GeV and V t bV t * = 0.04, we obtain the branching 
ratios as 


B(B S 


B(B S 


/iV) = 4.2 x 10“ 8 
= 6.8 x 10“ 9 
. T + T ~) = 8.9 x 10 -6 
= 1.4 x 10 “ 6 


(for 0 = 7 r) 
(for 0 = 0 ), 
(for 9 = 7 r) 
(for 0 = 0 ). 


(16) 


Thus, as seen from Eq. (HU) the branching ratio for B° s —*■ /r + /i has been enhanced by 
one order from its corresponding standard model value for 0 = 7 r, and is below the present 
experimental upper limit. This decay mode may be observable at the Tevatron Run II .is] 
to the level of 2 x 10 -8 . However, the predicted branching ratio for B® —> t + t~, which 
is C>(10 -6 ), could be observable in the currently running B factories, if we have a good 
efficiency for the detection of r lepton. 

Now let us consider the radiative dileptonic decay modes B s —> l + l~ 7 , which are also 
very sensitive to the existence of new physics beyond the SM. Due to the presence of the 
photon in the final state, these decay modes are free from helicity supression, but they are 
further suppressed by a factor of a. However, in spite of this a suppression, the radiative 
leptonic decays B s —> l + l~ 7 , l = (/r, r) have comparable decay rates to that of purely 
leptonic ones. The SM predictions for these branching ratios are B(B S —> /z + /i _ 7 , r + r~ 7 ) = 


1.9 x 10 , 9.54 x 10~ J respectively [17|,ll8|. These decays are also studied in some beyond 

the stanadard model scenarios jlf] |. 

The matrix element for the decay B s —> l + l ~7 can be obtained from that of the B s —> l + l~ 
one by attaching the photon line to any of the charged external fermion lines. In order to 
calculate the amplitude, when the photon is radiated from the initial fermions (structure 
dependent (SD) part), we need to evaluate the matrix elements of the quark currents present 
in between the emitted photon and the initial B s meson. These matrix elements can be 
obtained by considering the transition of a B s meson to a virtual photon with momentum 
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k. In this case the form factors depend on two variables, i.e., k 2 (the photon virtuality) and 
the square of momentum transfer q 2 = (p B ~ k) 2 . By imposing gauge invariance, one can 
obtain several relations among the form factors at k 2 = 0. These relations can be used to 
reduce the number of independent form factors for the transition of the B s meson to a real 
photon. Thus, the matrix elements for B s —> 7 transition, induced by vector, axial-vector, 


tensor and pseudotensor currents can be parametrized as 


20 ] 


{l(k,£)\s^^ 5 b\B s (p B )) = ie [e*(p B ■ k) - (e* ■ p B )k„] 


™ Bs 


('y(k,£)\sqf IJ ,b\B a {p B )) = ee^ a p£* u p% k 0 


Fy 

m Bs 


{'y(k,e)\sa IJU/ q v 'y 5 b\B B {p B )) = e • k) - (e* ■p B )k fl ] F TA , 

(^(k J £)\sa flu q u b\B s (p B )) = ee^ va0 £* v p a B k 0 F TV , 


(17) 


where £ and k are the polarization vector and the four-momentum of photon, p B is the 
momentum of initial B s meson and F t 's are the various form factors. 

Thus, the matrix element describing the SD part takes the form 

a 3 / 2 G F 


Msd = 


+ i(el(k-p B ) - (£* • p B )kn 


V tb V*\ e^ a ^p a B k 0 [A, Vfl + A 2 

£1 Ffi + B 2 h^ 5 l 


(18) 


where 


= 2 C^Ftv + C 9 , 
q 2 rn Bs 

B x = -2C 7 ^F ta - C 9 ^- , 
q m Bs 

The q 2 dependence of the form factors are given as 


Ao = C 


V 


10- 


™ Bs 


Bo = -C 


10- 


Fa 

m B 


F(E 7 ) = /3 


f Bs m Bs 

A + E’ 


where E 7 is the photon energy, which is related to the momentum transfer q 2 as 


E 1 = 


m Bs 


q z 


rrii 


(19) 


( 20 ) 


( 21 ) 


The values of the parameters are given in Table-1. 

When the photon is radiated from the outgoing lepton pairs, the internal bremsstrahlung 
(IB) part, the matrix element is given as 21] 


a 3 / 2 G F 

Mi B = — r— V tb V* f Bs mi Ci 0 

V 27 T 


J B 


J B 


p+ ■ k p_ ■ k 


75 l 


( 22 ) 
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TABLE I: The parameters for B —> 7 form factors. 


Parameter 

F v 

Ftv 

Fa 

Fta 

(3{ GeV" 1 ) 

0.28 

0.30 

0.26 

0.33 

A(GeV) 

0.04 

0.04 

0.30 

0.30 


Thus, the total matrix element for the B s —*■ l + l 7 process is given as 


M — M.sd + M-ib ■ 


(23) 


The differential decay width of the B —> / + 1 7 process, in the rest frame of B s meson is 


given as 


211 


dr g\ 


a 


ds 2 10 7t 4 


\VaV;,\‘ mi A , 


(24) 


where 


A = - s) 2y i((s + 2' r i)(\A 1 \ 2 + \Bi\ 2 ) + (s - 4ri)(\A 2 \ 2 + \B 2 \ 2 ) 

- 64 C\ 0 ((4 n - s 2 - 1) In | + V ' + 2s v t ) 

m Bs 1 — s V 1 — V[ J 

- 32 nil - sf f Ba Re(c 10 Aiy (25) 


with s = q 2 , s = s/m 2 Bs , 77 = m 2 /m 2 Bg , tq = \Jl — Arn 2 /q 2 . The physical region of s is 
4m , 2 < s < rn 2 Bg . 

The forward backward asymmetry is given as 


A fb 


1 

A 


2m 2 Bs s(l 


sfv 2 Re(AlB 2 + B*A 2 ) 


32 f Ba n{ 1 



Re(c 10 5 2 *) 


(26) 


Now using the form factors from m, we plot the dilepton mass spectrum (El, and the 
forward backward asymmetries (12611 for B s —> l + l ~7 decays which are shown in Figures-2 and 
3. In these plots we have used the weak phase difference between the SM and NP amplitudes 
6 to be 7 T to get the maximum possible contributions. From figures-2 and 3, we see that the 
branching ratio for B s —> l + l ~7 enhanced significantly from their corresponding SM values. 
However, the forward backward asymmetries are reduced slightly from the corresponding 
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FIG. 2: The differential branching ratio and the forward backward asymmetry ( Afb) for the 
process B s —> /r + /i _ 7 , in the standard model and in the NP model with the non-universal Z-boson 
effect. 




FIG. 3: Same as Figure-2, for the B s —> r + r 7 process. 

SM values and for the B s —> 7 process, there is a backward shifting of the zero position. 


To obtain the branching ratios it is necessary to eliminate the backgrounds, coming 
from the resonances with —» l + l~. We use the following veto windows to 
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eliminate these backgrounds 


B s -> /i + /i 7 : mj/^, - 0.02 < + 0.02; 

: m$> — 0.02 < m^+^- < m^> + 0.02 

B s —> r + r _ 7 : m^i — 0.02 < m T + r - < m^< + 0.02 . 

Furthermore, it should be noted that the \M.ib \ 2 has infrared singularity due to the emission 
of soft photon. Therefore, to obtain the branching ratio, we impose a cut on the photon 
energy, which will correspond to the experimental cut imposed on the minimum energy 
for the detectable photon. Requiring the photon energy to be larger than 25 MeV, i.e., 
Ery > 5 itibJZj which corresponds to s < m% s ( 1 — 6) and therefore, we set the cut 5 > 0.01. 

Thus, with the above defined veto windows and the infrared cutoff parameter, we obtain 
the brancing ratios as 

B(B° S -> fi+fi- 7 ) = 1.94 x 10 ” 8 , 

B(B° S -> t + t~ 7 ) = 1.37 x 10 ~ 7 , (27) 


which are enhanced by an order from their SM values. It should be mentioned that the 
B® —> t + t ~7 could be observable in the Run 11 of Tevatron. The contribution to the 
branching ratio due to bremsstrahlung photon is small for B s —> 7 and is found to be 

B(B° S —> /r + /i _ 7 )|ib ~ 0.5 x 10 -8 whereas it has dominant contribution to the B s —> t + t ~7 
process, i.e., B(B ° —»■ t + t~^)\jb ~ 1.3 x 10 - '. 

Now let us consider the flavor changing rare Z decays Z —> bs. Rare Z decays have been 
studied extensively in order to yield the signature of new physics. In the standard model 
this mode originates from one loop diagram with branching ratio ~ 3 x 10 8 . While the 

sensitivity of the measurement for the branching ratios for rare Z decays reached at LEP2 
is about 10 ~ 5 JrI, future linear colliders (NLC, TESLA) will bring this sensitivity up to 


10 -s level [23|. Various beyond the standard model scenarios has been employed [2J] where 
the branching ratio can be found to reach the sensitivity of the order of 10~ 6 . We would 
now be interested to analyze this decay mode in the model with an additional vector like 
down quark, where it can originate in the leading order. For completenes, we would also 
include the one loop corrections to the branching ratio, although their effect is negligibly 
small compared to the tree level contribution. 


10 




The tree level amplitude is given as 


M(Z —> bs ) = U sb 


9 


-U b ‘‘Yu Pl . 


2 cos 6w 

where denotes the polarization vector of Z boson. The decay width is is given as 


(28) 


r(z -> bs) = 


9 


rn z 


327r cos 2 6w 


\U sb \ 2 , 


(29) 


and the branching ratio to be 


B(Z —> bs + bs) = 


1 m z 


9 


-I U SI 


(30) 


T z 167T cos 2 6 W 

where T z is the total Z-boson decay width. 

Now we consider the one loop corrections arising from Figures 4-(b) and (c). The one-loop 
amplitude arising from Figure-4 (b) is given as 


1 


M{Z^bs) |i 

—loop 


167r 2 sb 2 cos 6^ 

+ B^rn^m^ml) + B 0 (m 2 ,m 2 ,m 2 z ) - C u (m 2 z , m 2 s , m 2 s , ml, m%, m 2 s ) 

+ 2C' 24 (m|, m 2 s , m 2 s , m 2 b , m|, m 2 ) + m| f 2Cn(m|, m 2 , m 2 , m 2 b , m|, m 2 ) 

+ 3C' 0 (m|, m 2 s , m 2 s , m 2 b , m 2 z , m 2 s ) - C 22 (m 2 z , m 2 s , m 2 s , m 2 b , m|, m 2 s ) 


9 2 

^b Pl 


O ri / 2 2 2 2 2 2\ 

2C 0 (m z , m s , m s ,m b ,m z , m s ) 


where 


, s-i ( 2 2 2 2 2 2 \ 

+ C 2 3 {m z , m s ,m s , m b ,m z , m s ) 


9 f 1 . 2fl 1 

flL = -T— q Sill - - 

COS 6V V d 4 


(31) 


(32) 


The contribution from Figure-4 (c) can be obtained from 4-(b) by replacing m s by m b and 
vice-versa. Including the one-loop correction the branching ratio for Z —> (bs + bs) is given 
as 

B(Z —■> bs + bs) = — —-—l^| 2 |l + -Ri + R\{jn s m;,)| 2 , (33) 

ix 167T cos- 





FIG. 4: Tree-level and one loop Feynman diagrams for Z —> 6s process. 
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where 


Ri = 


I67r 2 


s~i ( 2 2 2 2 2 2\ | n / 2 2 2 \ 

2 C'o(m z , m s , m s ,m b ,m z , m s ) + Bi(m s ,m s , m z ) 


iD/2 2 2 \ Py ( 2 2 2 2 2 2\. o/^~y / 2 2 2 2 2 2\ 

+ B 0 (m s ,m s ,m z ) - C n (m z , m s , m s , m b , m z , m s ) + 2C' 24 (m z , m a , m s , m b , m z , m s ) 
+ m| ( 2Cn(m|, m 2 s , m 2 s , m 2 b , m|, m 2 s ) + 3Co(m|, m 2 s , m 2 , m 2 b , m|, m 2 ) 


ry ( 2 2 2 2 2 2\ \ ry ( 2 2 2 2 2 2 \ 

- C 22 K, m s ,m s ,m b ,m z ,mj + C 2 z(rn z , m s ,m s ,m b ,m z ,mj 


(34) 


Using the quark masses (in GeV) as m 
the mass and width of Z boson from 
corrections as 



0.15 and m b = 4.4, sin 2 6 W = 0.23, a = 1/127, and 
, we obtain the branching ratio including one-loop 


B{Z -> bs) = 4.08 x 10 " 7 . 


(35) 


Since the expected sensitivity of giga -Z collider is of the order of 10 ~ 8 (which is at the level 
of SM expectation), we emphasize that new physics effect could be detectable in the rare 
decay Z —> bs, if indeed it affects this mode. 

Here, we have analyzed the rare decay modes B s —> l + l~ and B s —> l + l ~7 which are 
mediated by the b —> s FCNC transitions. We have considered the model which induces tree 
level FCNC coupling of Z boson, due to the addition of an extra vector like down quark to 
the matter sector. We found that the branching ratios for the radiave leptonic decay modes 
B s —> /i + /i~ 7 , (t + t“ 7 ) are of the order of 1CU 8 (1CU 7 ) which could be observable in the 
Tevatron Run II. Furthermore, the branching ratio of the pure leptonic mode B s —> t + t~ 
found to be 0( 10~ 6 ). This mode can be observed in the currently running B factories with 
improved r tagging efficiency. 

We have also analyzed the flavor changing decay of Z-boson to a pair of down quarks 
Z bs. This Z decay channel may prove useful in searching for new flavor physics beyond 
the SM at the TESLA or any other future collider which may be designed to run at the 
Z-pole with high luminosities, thus accumulating more than 10 9 on-shell Z bosons. With 
improved 6 -tagging efficiencies, the flavor changing decay Z —> 6 s is the most likely and the 
easiest one to detect among the flavor changing hadronic Z decays. It may be accessible to 
the Giga-Z option even for branching ratio as small as 13 (Z —> bs) ~ 1CU 7 — 10 -6 . 

To conclude, the standard model results of the rare decays studied here which are induced 
by FCNC transitions, are very small and cannot be detected in the current or near future 
experiments. These rare decays provide very sensitive probe of new physics beyond the SM. 
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Detection of these decays at visible levels by any of the future colliders would be a clear 
evidence of new physics. 
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APPENDIX A: ONE-LOOP FORM FACTORS 


The two-point and three-point one-loop form factors which are defined as 


C 0 ; C 

_ f d 4 q 1; qpq» 


m 1 2 3 [q 2 — mf] [(q + pi) 2 — m|] [(q — p^) 2 — m|] 


(Al) 


Co-.C^lm^m^m^p^p^p^) = 


d A q 


q 2 ; 


m 2 [q 2 — ml] [(q + pi) 2 — ml] [(q — ps) 2 — ml] 
where "YhiPi = 0 is to be understood above. 


, (A2) 


Bo; B^im^ m 2 ,p ) = 


d 4 q 


1 ; Qa 


in 2 [q 2 — ml] [(g + p) 2 — m 


21 ’ 


(A3) 


The coefficients Bj with j G 0,1, Cj with j G 0,11,12, 21, 22, 23, 24 are defined through 
the following relation 


B IX = p fl B 1 

Cfi = PlfiCli + P2fj, + P2fiCl2 

Cfiu = PliiPlvC'2i + P2fiP2uC22 + {Pi P2 }ilvC23 + SVt/CW 

C, iy = Pl^Pl v C21 + P2 f ,P2uC 2 2 + {Pi P2}^C 2 3 + P^C 2 4 (A4) 

where {a b} = + a u b A1 . 
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